NCODA  Status 

NRL  Coupled  Ocean  Data  Assimilation 

James  A.  Cummings 
Naval  Research  Laboratory 

HYCOM  NOPP  GODAE  Meeting 
RSMAS,  University  Miami,  6-8  December  2005 

Outline 

1.  NCODA  System  Overview 

2.  New  Analysis  Capabilities 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
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NCODA  System  Overview 


Flexible  System 

-  global  or  regional  applications 

-  re-locatable,  multi-scale  analyses  on  nested,  successively  higher 
resolution  grids  (3:1  nest  ratios) 

-  update  ocean  forecast  model  or  run  stand-alone 

•  2D  analyses  of  sea  ice  and  SST  (NWP  boundary  conditions) 

•  3D  temperature  and  salinity  analysis  (geostrophic  currents) 

•  3D  MVOI  sequential  incremental  update  cycle  (model-based) 

Designed  as  Complete  End-to-End  Analysis  System 

-  data  quality  control,  analysis,  performance  diagnostics 

-  operational  at  Navy  Centers  in  analysis-only  mode 

•  Naval  Oceanographic  Office 

•  Fleet  Numerical  Meteorology  and  Oceanography  Center 


Ocean  Obs 


NCODA  Implementation 

Sequential  Incremental  Update  Cycle 
Analysis-Forecast-Analysis 


SST:  Ship,  Buoy, 

AVHRR  (GAC/LAC), 

GOES,  AMSR-E, 

MSG,  AATSR 
SSS:  TSG 
Temp/Salt  Profiles: 

XBT,  CTD,  Argo  Float, 

Buoy  (Fixed/Drifting), 

Gliders 

SSH:  Altimeter,  T/S 
profiles 

Sea  Ice:  SSM/I  Forecast  Fields 

Prediction  Errors 


Innovations 


3D  MVOI 


Increments 


First  Guess' 


Ocean  Model 


Model  forecast  fields  and  prediction  errors  are  used 
in  the  QC  of  newly  received  ocean  observations 


New  NCODA  Capabilities 


Analysis 

-  first  guess  appropriate  time  (FGAT) 

-  flow  dependent  correlations 

-  analysis  error 

-  model  climate  error  variance  fields 

-  age  of  data  on  grid  analysis  variable 

-  pressure  correction  analysis  variable 

-  data  restriction  in  boundary  areas 

Observing  Systems 

-  MeteoSat  SST  (MSG)  from  MeteoFrance  (CMS-Lannion) 

-  Microwave  SST  (AMSR-E)  from  Remote  Sensing  Systems 

-  Ocean  Gliders  (up/down  profiles,  position  varies  with  depth) 


AVHRR-GAC  AMSR-E 
12-hrs  data 


Sequential  Incremental  Update  Cycle 
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Innovations 


Length  update  cycle  user  defined 

All  observations  considered  synoptic  regardless  length  update  cycle 


First  Guess  at  Appropriate  Time 
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Data  Window  (+/- 12  hours) 
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Innovations 


Length  update  cycle  user  defined 
Interval  of  forecast  periods  user  defined 
Eliminates  component  of  mean  analysis  error  that  occurs  when  comparing 
observations  and  forecasts  not  valid  at  same  time 


Flow  Dependent 
Correlations 

hs  =  0.2 

small  (large)  hs  produces 
strong  (weak)  flow 
dependence 


Sh=(x0-xb)/xs  C„  =(l  +  sfc)exp(-sj 

sv=(z„-zb)/ zs  Cv  =(l  +  sv)exp(-sv) 

sf  =(h0-hb)/hs  Cf  =  (1  +  sf )  exp(-s/ ) 

c  =  c  c  c 

xs  =  horizontal  scale  (km)  Ch  =  horizontal  correlation 
zs  =  vertical  scale  (m)  Cv  =  vertical  correlation 

hs  =  flow  scale  (dyn.  m)  Cf  =  flow  correlation 


Geopotential  Valid  5  August  2005  00Z 

Contour  Interval  0.1  dyn  m 
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100  M  Temperature  Increments 
6  August  00Z 


Analysis  Error  Reduction  (%) 


Satellite  SST  + 

In  Situ  SST 
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GAC  D/N  GOES  D/N  LAC  D/N 


2fc  Ship  E\red  Hucit  Drft  Bu ay  Sea  Ice  33V  Climate  CKAK  3J7I  toaljaed  33T  .aircraft  33T 


Pa=Pb-PbHT[HPbHT+R]'HPb 

estimates  the  forecast  error  reduction 
due  to  the  observing  networks 

Pa  -  analysis  error 
Pb  -  background  error 
R  -  observation  error 
H  -  measurement  functional 


0  20  40  60  80  100 

Error  Reduction  (%) 

Will  be  used  in  ETKF  to  determine 
impact  of  glider  data  assimilation 
(adaptive  sampling  RTP) 


Model  Based  Error  Variances 


•  computed  from  differences  of  free  running  model  states  at  analysis 
update  cycle 

•  provides  estimates  of  model  error  (variability)  for  all  analysis  variables 
(T,S,u,v,h) 

•  used  by  NCODA  in  time  evolution  of  background  error  variances 


C^™04'  '  .OB '  '  .'is'  ' 

HYCOM  SSH  Model  Climate 
Variability  -  9  km  grid  Gulf  Mexico 


NCODA  Background  Error 
Variances 

•  vary  by  position,  depth,  analysis  variable 

•  evolve  with  time,  updated  continuously 
using  analyzed  increment  fields 

•  error  growth  parameterization  in  data 
void  areas 

-  function  of  age  of  data  on  grid  and 
temporal  autocorrelations 

-  background  errors  asymptote  at 
model  (climate)  variability  in  long 
term  absence  of  observations 


New  Analysis  Variable:  Age  of  Data  on  Grid  (hrs) 

number  hours  since  grid  point  influenced  by  an  observation 


In  Situ  SST 


Age  of  Data  at  Surface 


MODAS 

Synthetics 


Age  of  Data  at  400  M  Depth 
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SSH  Increment  (m) 
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Pressure  Correction  (db)  at  100  M 


New  Analysis  Variable: 
Pressure  Correction  (db) 

compute  pressure  innovation  (OmF)  of 
forecast  density  in  observed  density  profile 

compute  pressure  innovation  error  from 
T,S,p  errors  scaled  by  observed  potential 
density  pressure  gradient 

ep  =(e0  •  dp/d6  +  es  •  5p/<9S)  /  (dp/ dp) 

assimilate  pressure  innovations  using  T/S 
covariances  in  T,S,(|)5u5v5p  analysis 

correct  HYCOM  forecast  interface  pressures 
when  layers  are  at  target  density 

-  positive:  move  the  layer  down,  forecast 
density  shallower  than  observed 

-  negative:  move  the  layer  up,  forecast 
density  deeper  than  observed 

correct  HYCOM  forecast  T,S,p  when  layers 
are  not  at  target  density 

apply  constraints  before  initializing  model 

-  layer  thickness  is  always  positve 

-  no  bottom  pressure  change 


Data  Restriction  in  Lateral  Boundary  Areas 


Topex  ERS2  GFO  Jason  Envisat  In  situ  Suppl 

Altimeter 
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Satellite  SST 
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Argo  Float  FLsed  Buoy  Drlltlng  Buoy 


In  situ  SST 


Profiles 


•  model  forecast  not  accurate 
in  lateral  boundary  areas 

•  innovations  (OmF)  can  be 
large  in  boundary  areas 

•  analyzed  increments  from 
boundary  areas  can  degrade 
forecast  in  active  regions 


Temperature  Increments 
200  M 
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